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1. Introduction 



Suppose that a mass standard having a small but finite 
magnetic susceptibility or permanent magnetization has 
been placed on one pan of a beam balance. The mag- 
netic properties of the standard may lead to a potential 
energy term that also depends on the magnetic fields 
within the balance. If this magnetic potential energy 
changes as a function of rotation of the balance beam, 
then there will be an unwanted torque on the beam. One 
may also note that rotation of the beam through small 
angles is equivalent to motion of the mass standard in the 
vertical direction. From this qualitative argument, it is 
clear that a magnetic potential due to the mass standard 
and any balance components fixed with respect to the 
frame of reference of the standard can pose no problems 
to weighing. Problems, when they occur, are due to the 
magnetic properties of the mass standard and local mag- 
netic fields that are in a frame of reference (usually that 
of the laboratory) that moves relative to the standard on 
the balance pan. In order to quantify such problems it is 
necessary to characterize both the magnetic properties 
of the mass standard and the magnetic environment of 
the laboratory. This paper will concentrate on the 
former task. 



In a previous publication [1] we showed how the 
volume magnetic susceptibility x of nonmagnetic mate- 
rials could be measured by simple modification of a 
commercial microbalance. Tests carried out at the 
Bureau International des Poids et Mesures (BIPM) veri- 
fied that a reasonably uncomplicated theory of opera- 
tion is adequate to achieve results that have a relative 
combined standard uncertainty of a few percent over the 
range of susceptibilities encountered in mass standards 
of high quality. (Throughout this report uncertainties are 
expressed according to guidelines set forth by the Inter- 
national Organization for Standardization [2].) 

In this report, we first discuss a model for errors in 
mass metrology due to magnetic effects. These are 
shown to depend mainly on the susceptibility and the 
permanent magnetization of the mass standard. The 
latter is zero for nonmagnetic materials but some mate- 
rials used in the manufacture of good mass standards 
are, in fact, weakly magnetic and thus may become 
magnetized. 

Susceptibility can be measured quantitatively using 
the BIPM apparatus and permanent magnetization can 
be detected. These measurements are discussed in 
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detail, including calibration of the instrament and esti- 
mation of measurement uncertainty. Of equal impor- 
tance to making good measurements is interpreting their 
significance. We therefore begin with a basic review of 
the most probable magnetic problems in precision mass 
measurement and conclude with several practical exam- 
ples. 

2. Model for Magnetic Errors in Mass- 
Metrology 

By magnetic error, we mean an unsuspected vertical 
force F that is magnetic in origin. Such a force will be 
misinterpreted as a mass Fig , where g is the local accel- 
eration of gravity. 

We may assume that high-quality mass standards are 
artifacts with an isotropic volume magnetic susceptibil- 
ity X whose magnitude is much less than one. In addi- 
tion, artifact mass standards should have little or, ideally, 
no permanent magnetization. Thus the unwanted mag- 
netic force will, to a good approximation, be given by 



t!±A. 

' 2 dz 



jx'H-HdV-f,o-l-J 



MHdV, 



(1) 



where x' is the effective volume magnetic susceptibility 
of the standard, M is its permanent magnetization (de- 
fined as the magnetic moment per unit volume in zero 
field), H is the local magnetic field strength and the 
z-axis is parallel to g. In general, all these terms may 
depend on position but we assume that ;if ' is a scalar. The 
parameter /Xo is the vacuum permeability, identically 
equal to 4ttX 10'^ N- A"^. The integrals are taken over the 
volume of the artifact. The effective susceptibility x' is 
defined as X~Xa, where Xa is the volume susceptibility 
of air (+3.6X 10''). Because of its relatively small mag- 
nitude, Xa can be neglected in most of the examples 
given below. The symbols, quantities and nomenclature 
used in this report are those set forth by the Interna- 
tional Union of Pure and Applied Physics [3]. 

We have assumed in Eq. (1) that H is the field before 
the sample is introduced, an approximation good to first 
order in the susceptibility. We have also assumed that 
the alloy is linear, i.e., its susceptibility is independent of 
applied magnetic field for strengths less than, say, 
4 kA-m"'. Finally, we have assumed that the effect of a 
permanent magnetization M can simply be added as a 
term separate from the induced magnetization x'H- The 
validity of this approach will be demonstrated by exam- 
ples given below. Note that, for a mass standard to be 
magnetized, it must have been exposed to field strengths 
greater than the linear threshold. 

We can further simplify Eq. (1) by the additional 
assumption, usually warranted, that the linear and 



isotropic susceptibility is also homogeneous throughout 
the artifact. Finally, we make an mm warranted assump- 
tion that M is a constant in magnitude and direction 
throughout the artifact and is independent of H at low 
field strength. It would be difficult to proceed without 
the last assumption but the reader should keep in mind 
that its importance in what follows is largely heuristic. 
This allows us to write a simplified force equation: 



F = 



IM> 



2 ^ dz 



HdV - iXoM, 



fJioM, 



dz 






/. 



dV 



(2) 



It should be noted that M^, M^ and My are the vertical 
and horizontal components of M and thus change with 
the orientation of the artifact. Without knowledge of the 
magnetic fields and gradients within the balance our 
model can take us no further except to imply the follow- 
ing: (i) There can be no magnetic errors if there are no 
field gradients in the balance, and (ii) Due to the sym- 
metry of certain shapes (e.g., cylinders) some artifacts 
can be turned upside-down without changing the force 
contribution from the susceptibility. The sign of M^ is, 
by contrast, reversed for these shapes. 

Two serious attempts have been made to apply Eq. (2) 
to problems of mass metrology. Gould [4] presented 
several special cases to illustrate his contention that the 
worst problems are likely to occur when the mass stan- 
dard and a part of the balance case (including the mass 
transporter in automated balances) located just below 
the pan are both accidentally magnetized along the ver- 
tical axis. This argument implies that the second integral 
in Eq. (2) is the most important contribution, the fields 
and gradients arising primarily from the magnetized 
parts of the balance. 

Using this reasoning as a guide for the selection of 
suitable stainless steel alloys, Gould concluded that an 
alloy used for mass standards should be chosen both for 
its low susceptibility and for its resistance to permanent 
magnetization upon exposure to high fields. He found 
that for stainless steels, these two desirable properties 
are well correlated. That is, the alloy with the lowest 
susceptibility was also the most difficult to magnetize. 
The recommended alloy was found to have x ~ 0.003 
when measured in a uniform field strength of 
16 kA-m"'. The permanent magnetization was less than 
1 A-m"' after exposure to a "suitably large" uniform 
field (we have converted values given in Ref. [4] from 
CGS-EMU to SI [3]). 

Kochsiek [5] approached the problem in a different 
way. He dealt with Eq. (2) by arguing that, once demag- 
netized, the normal use of mass standards should not 
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subject them to fields great enough to remagnetize 
them. A strong recommendation was made against se- 
lecting inferior stainless steel alloys, known to be easily 
magnetized. He then chose a special case in which the 
integrals in Eq. (2) may be solved trivially: the magnetic 
field within the balance chamber is vertical with the 
form 

H, = ho + hiz, (3) 

where ho and hi are constants and the origin of the 
coordinate system is at the geometric center of the mass 
standard. As an example, ho was chosen to be 100 A-m"^ 
and hi to be 5000 A-m"^. (As a point of reference, the 
vertical component of the Earth's magnetic field 
strength is about 40 A-m"' at the latitude of Paris.) 
Based on this somewhat arbitrary choice of "worst- 
case" parameters, it was then possible to argue that 
secondary 1 kg mass standards used by national labora- 
tories should have a volume susceptibility below 0.003. 

It is instructive to note that this recommendation also 
allows us to calculate an upper limit for M^ in Eq. (2): 
IMJ < 0.3 A-m"'. The limit is derived from Eqs. (2) and 
(3) with the condition that the second integral in Eq. (2) 
be smaller than the first when x' = 0.003. In general for 
this model, M^ is negligible when \M^ I < I x'ho I- 

Thus we see that Refs. [4] and [5], while focusing on 
different aspects of the problem, arrive at nearly the 
same guidelines for selecting stainless steel alloys suit- 
able for the highest quality of secondary mass standards. 
Nevertheless, it should be emphasized that Eq. (2) can- 
not be solved without detailed knowledge of the fields 
within the balance. The Organisation Internationale de 
Metrologie Legale (OIML) has recently recommended 
that Class Ei and E2 mass standards have volume suscep- 
tibilities less than 0.01 and 0.03 respectively [6]. 



3. Susceptometer Developed at the BIPM 

The BIPM susceptometer subjects a mass standard to 
relatively small, calculable fields, measures the resulting 
vertical force and then makes use of Eq. (2) to find the 
susceptibility and the parameter M^ . The principles of 
construction are given in Ref. [1] but we show the ap- 
paratus schematically in Fig. 1 and give design details in 
Appendix A. The magnet that we use is cylindrical with 
height L and diameter both equal to 5 mm. 

The balance reads in units of mass, the reading be- 
coming more negative when a paramagnetic sample is 
introduced. Thus, the following relation is used : 

-Hi (xN corresponds to a balance reading of -101.9 (xg, 

where the acceleration of gravity g has been taken to be 
9.81 m-s"'. 



©. 



H© 



®M 



&^ 



^><S) 



Fig. 1. Schematic view of the apparatus. A small rare-earth magnet 
(A) is placed on a supporting column which rests on the balance pan 
(B). A nonmagnetic bridge (C), the height of which may be incre- 
mented with gauge blocks (D), straddles the balance. The sample (E) 
is placed on the bridge, directly above the magnet. 



3.1 Response to Large Samples 

Some insight can be gained by considering a sample 
in the form of a semi-infinite slab and approximating 
the magnet as being uniformly magnetized along its 
axis. If the magnet were a sphere, the first integral in Eq. 
(2) would have a particularly simple solution [1]: 



F = 



3fjioX'm^ 
64 ttZS 



■^X'Fr. 



(4) 



where m is the moment of the magnet (Sec. 4) and Zo is 
the distance between the center of the magnet and the 
sample, as shown in Fig. 1. The remaining integrals 
vanish in this case. A more general result, correct to all 
orders of x, could have been derived using the method 
of images. 

The method of images may also be applied to the case 
of the cylindrical magnet actually used, modeling the 
magnet as a cylindrical sheet of uniform current. The 
solution involves elliptic functions that can be calculated 
numerically. Only the first integral of Eq. (2) is not zero, 
just as in the case of a spherical magnet. In Fig. 2, force 
calculations for cylindrical magnets of different aspect 
ratio y (height/diameter) are compared with calculations 
for a spherical dipole, where all magnets have the same 
volume and the same value of uniform magnetization. 
The sample is again assumed to be a semi-infinite slab 
of small susceptibility. It is, of course, well known that 
cylindrical magnets cannot be uniformly magnetized 
but the approximation is a good one for the neodymium- 
iron-boron magnets which we use. 

For large samples, a cylindrical magnet with y = 1 
thus produces almost the same force as a spherical mag- 
net of equal moment. The difference is within 2.6 % for 
Zo>2L and within 1.2% for Zo>3L. Even closer 
agreement would be expected if y were 0.87 instead of 
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1. However, the usual shape of cylindrical rare-earth 
magnets has y = 0.5, so it is convenient simply to com- 
bine two such magnets. 

Using the property of superposition, we can immedi- 
ately find the force on a slab of finite thickness i . If we 
define Zi = Zq + i, then 



F = X'Fm 



1 



X'Frn 



1-Zl 



where the dimensionless quantity Zi„ is the length Zi 
normalized to Zq. 

More generally, we next consider how large a cylin- 
drical sample of radius c and height t must be before it 
may be approximated as "semi-infinite." To do this, we 
insert the dipole fields due to the magnet in the first 
integral of Eq. (2) and solve for the case of cylindrical 
samples of finite dimensions. The result of the calcula- 
tion is given in Eqs. (5) and (6). 

Figure 3 shows contours of equal force between the 
magnet and samples of differing dimensions. The force 
is normalized to x'Fmax and the cylinder dimensions are 
normalized to Zq. Only a finite volume of semi-infinite 
samples contributes strongly to the force integral. The 
important volume is of order Zo^ and thus changes as Zo 
is varied experimentally. 
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Fig. 2. Expected behavior of three cylindrical magnets of different 
aspect ratio (height/diameter) y relative to a spherical magnet. All four 
magnets have the same volume (98 mm') and uniform magnetization. 
The case considered is the attractive force between each magnet and 
a semi-infinite slab of small susceptibility. The spherical magnet is 
taken as the reference and is thus represented by the abscissa. 
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Fig. 3. Contours of constant force between a dipole magnet and coaxial cylinders of 
varying radius c„ and thickness f„ but having the same, small susceptibility. Cylinder 
dimensions are normalized relative to Zo, the distance between the center of the magnet and 
the base of the cylinders. Force contours are normalized to the signal from a cylinder of 
infinite radius and height (i.e., a semi-infinite slab). 
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3.2 General Equations 

A sample placed on the susceptometer as shown in 
Fig. 1 will be subjected to fields from the cylindrical 
permanent magnet (//mag) and the Earth (//e)'- Making 
the dipole approximation for the cylindrical magnet and 
taking the Earth's field to be uniform leads to the fol- 
lowing force equation [1] 



F = x' ^max h + j^ (x'He, + M,) ^ /b = F, + Fb , (5a) 

477 Zo 

where 



32 
3m 



? AJ ]]"-''' "-'''^^'- 



ti-Jii^^^'^'''^^ '''' 



slightly each time the magnet is repositioned. In prac- 
tice, we have not found this to be a problem for nominal 
settings of about 15 mm or greater. 

From a measurement of F^ using the unknown sample 
and a knowledge of m, Zo and the sample dimensions, 
we can solve for x\ using Eqs. (4), (5a) and (5b). Using 
this x\ the measured value of Fy, , and a knowledge of 
Hez we can then find M^ from Eqs. (5a) and (5c). The 
uniform magnetization is usually negligible. The quanti- 
ties /a and /b are readily written in closed form for 
cylindrical samples of radius c and height t coaxial with 
the magnet. This is because each differentiated triple 
integral reduces to a single integral in this case. We 
repeat the relations, given in Ref. [1], using normalized 
variables^: 



/a=l-Zi; 



1 +cl/3 



+ Zv 



1 + (cJZnf/3 

[1 + (cJZ,f]' 



(6a) 



and 



477 _a_ 

m dZo, 



ill 



(/fmag)z dV-- 



277 



Cn/Zln 



(i+ciy^ [i + (c„/Zi„)T 



(6b) 



d 

dZon 



/// 



■2z' 



ip' + z^' 



pApAOAz 



(5c) 



and Fmax is defined in Eq. (4). We have assumed the 
samples to be centered symmetrically about the magnet 
axis so that only the vertical components of //e and M 
contribute to the force. The integrals are taken over the 
normalized volume of the sample. It is evident from Eq. 
(5a) that the induced magnetization x'Hez has the same 
effect as a permanent magnetization M^ . We show below 
how the two contributions may, in practice, be distin- 
guished. The differentiated integrals of Eqs. (5b) and 
(5c) are shown in detail for cylindrical coordinates 
where the origin is, as always, at the center of the mag- 
net. 

In practice, an initial force measurement Fi is made 
with the north pole of the magnet pointing down and a 
second measurement F2 is made (at the same Zo) with 
the north pole pointing up. Then the first term in Eq. 
(5a) is given by F^ = (Fi + F2)/2 and the second term by 
Fb = (Fi - F2)/2. It was noted in Ref. [1] that Zo changes 



where Zi„ = 1+ t„. Note that < I^< 1 , the upper limit 
being the case of a semi-infinite slab and the lower limit 
a sample of vanishingly small volume. In contrast, 
s /b s 2.42, the lower limit being obtained for both a 
semi-infinite slab and a sample of vanishingly small 
volume. For shapes other than cylinders, solutions can 
be found by numerical integration and/or superposition 
of easily calculable shapes. Several examples are given 
below and in Appendices B and C. 
It can be shown in general that 

Fa _ 377 X'timax h_ 

Fb ~ 8 x'He. + Mz h ' 



where 



tlm 



IttZV 



the magnitude of the maximum field strength to which 
the sample is subjected by the magnet alone. This is the 
field at the base of the sample, directly above the mag- 
net. An interesting special case, which is not unusual, is 



' This assumes that the laboratory is not magnetically shielded and 
that no extraneous fields are present. These assumptions were verified 
for our laboratory. 



" Equation (8) of [1] has obvious typesetting errors. 
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that of a cylinder for which Z^^ » 1 and iZ^Jc^ » 1. 
Then 4 takes its maximum value at c^ = v2, in which 
case /a//b = 0.4. Thus we see that, for a given Zo, an 
apparatus realized with a much more sensitive balance 
and a much smaller magnet than we use would favor the 
measurement of F,, instead of F^. 

At any given setting of Zo, the geometric terms in Eq. 
(5) are the same for samples having the same dimen- 
sions and orientation. Thus the relative susceptibilities 
of congruent samples are easily determined from ratios 
of balance readings. To determine an unknown suscepti- 
bility, one might even machine a material of known 
susceptibility to match the dimensions of the unknown 
object [7]. Unless a large number of identical samples is 
to be determined, or the field of the magnet cannot be 
approximated as a dipole, this approach is less conve- 
nient than finding sufficiently accurate solutions to the 
differentiated integrals of Eqs. (5b) and (5c). 

The susceptibilities of a set of objects, each having 
the same shape and size (such as similar mass stan- 
dards), may be conveniently determined by determining 
the susceptibility of one and then treating it as a stan- 
dard to find the susceptibility of the others through 
ratios of the observed values of F^. 

It now remains to show how the BIPM susceptometer 
may be calibrated and how we use it. 

4. Calibration of Susceptometer 

Calibration of the device is equivalent to determining 
m and Zq. Several methods may be used and they give 
equivalent results. The two unknown parameters may, 
for example, be found by measuring two standards of 
different, known susceptibility. Such standards are not 
readily available in a suitable range of susceptibilities, 
however [8]. A reasonably good estimate of m can be 
inferred from technical information supplied with the 
magnet. As shown below, the susceptometer may also 
be self-calibrated by fitting measured values of F^ as a 
function of increments in Zo , obtained using a sample of 
unknown susceptibility, to the equations developed 
above. All these methods should, of course, be consis- 
tent within their combined uncertainties. 

4.1 Determination of m From Supplier's Specifica- 
tions 

At distances Zq that are large compared with the mag- 
net dimensions, the axial magnetic field strength due to 
the magnet is given by //^ax, and may be positive or 
negative depending on the orientation of the magnet. If 
the magnet has a uniform axial magnetization M, then 
m = MY where Y is the volume of the magnet. It is more 



usual to characterize a permanent magnet by its polar- 
ization / = /UoM, which has units of tesla. 

The supplier of our magnets provided a depolariza- 
tion curve that shows how / is reduced in the presence 
of an opposing external field (Fig. 4). The shape of the 
magnet also has an effect on /. Even with no applied 
field, a cylindrical magnet will suffer to some extent 
from self-demagnetization if y < oo. In the case of a 
cylindrical magnet with 7=1 (height equal to diame- 
ter), the demagnetizing factor A^^ is 0.312 [9]. Thus a 
first estimate of m is the volume of the magnet multi- 
plied by the value of M corresponding to the dot in Fig. 
4. Note that the demagnetizing factor of the magnet has 
a relatively small effect on M. This is an important 
property of rare-earth magnets. 




J/T 



Fig. 4. Depolarization curve supplied by the manufacturer of the 
magnet used in the susceptometer. The ordinate represents the mag- 
netic polarization and the abscissa the magnetic field strength oppos- 
ing the polarization. For a cylindrical magnet with height equal to 
diameter (A'„ = 0.312 [9]), self-demagnetization shifts the operating 
point from the /-intercept (no opposing field strength) to the point 
indicated by a dot. 



Taking M and Y from the manufacturer's data sheets 
thus gives us the following estimate: m = 0.0898 A-m^. 
It is difficult to assign an uncertainty but we will see in 
Sec. 4.2.2 that experimental values are within 1.5 % of 
this. 

The estimated value of m is useful in computing an 
approximate value for //^ax, the maximum field strength 
to which the sample will be exposed at a nominal setting 
of Zo. For 15 mm <Zo<30 mm, 4.25 kA-m"' > 



//„ax > 0.53 kA-m-\ 
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4.2 Bootstrap Calibration 

The calculation of m given in Sec. 4.1 is a good 
approximation but the definitive value for m is best de- 
termined experimentally. In any case, we still require a 
routine method of finding Zo to sufficient accuracy. Both 
needs can be met by carrying out an initial series of 
bootstrap measurements, shown schematically in Fig. 5 
and described in detail below. We first determine x'^^^ 
for a selected standard S. Based on Eqs. (4) and (5a), this 
quantity then allows us to find Zo by placing the standard 
on the apparatus and measuring F^. Once we have a way 
of determining Zo, we can find m, as shown below. The 
bootstrap measurements require considerable effort but 
need only be done once. In the following, we analyze 
statistical uncertainties (i.e.. Type A standard uncertain- 
ties) only. Type B standard uncertainties are discussed in 
Sec. 4.3. 

4.2.1 Determination of x'sin^ (and hence 
Zo) Using gauge blocks (see Appendix A), Zo can be 
incremented in precise steps of, for example, 5 mm. The 
value of Zo may thus be known to within a constant 
length Zoo = Zo - Zb, where Zb is the height of the gauge 
blocks in use. If we had a suitable standard with known 
susceptibility x^ ^nd if, further, we knew m^, then a 
measurement of F^ using the standard would determine 
the value of Zo each time the height of the stage was 
changed. 



Standard samples of sufficiently large susceptibility 
are notoriously difficult to obtain [8] and so we stan- 
dardize an unknown sample by measuring it at different 
settings of Zo and fitting the results to the function 
predicted by the theory presented above. The parame- 
ters obtained from this fit contain some of the informa- 
tion we require. The force F^ computed from Eqs. (4), 
(5a) and (6a) depends on the dimensions of the standard 
as well as the quantities (x'stn^) and Z). For a suitable 
standard which is linear (up to field strengths of about 
20 kA-m"'), isotropic and homogeneous, F^ as a func- 
tion of Zb may be fitted using two adjustable parameters 
and hence the unknowns x'^fn^ and Zo can be found. (In 
the following we will use the abbreviation LIH to refer 
to materials that are magnetically linear, isotropic and 
homogeneous.) As long as we use the same magnet and 
the same sample, ;^'sm^ will not, in the normal course of 
events, change with time and so the standard may be 
used to find Zo whenever the height of the apparatus is 
changed. 

Given a cylindrical sample that is LIH and has a 
susceptibility of about 0.0015, the major difficulties are: 
i) The computation of /a requires a value for Zoo, which 
is unknown. Thus all calculations must be carried out 
iteratively by making an initial guess of Zoo. 
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Fig. 5. Flow charts showing the two steps of self-calibration. Step 1 determines the 
quantity Xa"^^ ^nd Step 2 then determines the individual parameters m and x's- 
Computations shown in Step 1 can be greatly simplified through the use of commer- 
cial software. 
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ii) The functional form of F^ versus Zb is not linear. 
A linear result will be obtained by plotting (h/F^y^ 
versus Zb. 

iii) As shown in Sec. 3.1, results obtained at relatively 
small values of Zo are biased by assuming the cylindri- 
cal magnet to be a dipole. 

The experimental algorithm is shown as Step 1 of Fig. 
5. The sample we chose is made of the nonmagnetic 
alloy Alacrite X.S.H. (Aubert and Duval S. A., Neuilly 
sur Seine, France)^ whose nominal mass composition is: 
20 % Cr, 15 % W, 10 % Ni, 0.1 % C, remainder Co. It 
is a polished disc with a diameter of 69 mm and a 
thickness of 9.8 mm and had already been fabricated for 
another purpose. In retrospect, a thickness of about 
20 mm would have been preferable to increase the sig- 
nal. Measurements of F^ were made starting at Zq = Zoo 
~ 9.5 mm and in precise increments of 5 mm up to Zq = 
Zoo + 20 mm. 

Table 1 shows the final iterated values of ;^''sm^ and 
Zoo based on three different least squares fitting routines. 
For fit (1), a weighted linear least-squares routine was 
used, where the uncertainty in Zb is assumed to be 
negligible. Formally, we are fitting y = l3\-x+l3o where 

/ / \ 1/4 

J =1 — 1 and x = Zb. Since /^ depends weakly on 
iPJ Pi)+x, several iterations are required. 

Table 1. Data and data analyses for determination of the quantity 
)('sm^ from measurements of F^ as a function of Zb. Standard uncer- 
tainties for the measurements of f „ are given in the second column 



f./ynN" 




Mc//xN 




Zfl/mm 


24.65 




0.07 







3.990 




0.0012 




5 


1.112 




0.010 




10 


0.4042 




0.005 




15 


0.1702 




0.005 




20 


Fit number'' 




1 




2 


3 


X\m^l{KW) 10.718X10- 


"■ 10.539X10- 


"■ 10.703X10^ 


o-/(AW) 0.069X10- 


' 0.10X10-' 


0.044X10-' 


Zoo/mm 


9.477 




9.404 


9.470 


o-/mm 


0.024 




0.075 


0.014 


Deg. freedom 


2 




2 


2 



" +1 jjlN corresponds to a balance reading of -101.9 (xg. 
"" The three fits shown in the table result from different ways of manip- 
ulating the same data. Fit (3), which was carried out using commercial 
software, is taken as the definitive result. 



Certain commercial equipment, instruments, or materials are identi- 
fied in this paper to foster understanding. Such identification does not 
imply recommendation or endorsement by the National Institute of 
Standards and Technology, nor does it imply that the materials or 
equipment identified are necessarily the best available for the purpose. 



The datum at Zo ~ 9.5 mm, being the most biased by 
use of the dipole model and also being outside our 
normal range of use, was simply not used in the fit. Fit 
(2) differs from (1) in that equal weight was given to the 
four estimates of iUF^)"'*. Use was made of the covari- 
ance matrices for each fit in order to arrive at the tabu- 
lated standard uncertainties. We expect the equal 
weighting of fit (2) to yield unbiased estimates of the 
regression parameters but with a larger uncertainty 
compared to a weighted fit. The tabulated values con- 
firm this expectation. 

Fit (3) uses the curve-fitting utility of Sigmaplot 1 .02 
to fit the data directly to Eqs. (5a) and (6a). Formally, 

we are fitting y = — — '-^-^ where j = Fa , x = Zb , and I^ 
(Ao + x) 

is an implicit function of Ao -n x . The results shown in 
Table 1 are for equal weight given to the data. Results 
obtained by weighting the input data are scarcely differ- 
ent. Again the point at Zb = has been omitted. The two 
parameters shown in the table are highly correlated in all 
three fits. 

Taking fit (3) as our estimator, the final result is 

X's m^ = (10.70 ± 0.05) X 10-'^ A'•m^ 

We note that using the fitted values to extrapolate to Zb 
= results in a value of F^ that is about 5 % lower than 
the observed force. The implications of this error are 
discussed in Sec. 4.3. 

4.2.2 Determination of m Given ;^ ^ m , it is a rel- 
atively simple step to determine m , and hence x's (Step 
2 of Fig. 5). To do this, we need two additional magnets, 
B and C, of similar dimensions to A, the magnet we 
normally use. The additional magnets have moments 
that are, as yet, undetermined although all three mo- 
ments are expected to be consistent with the estimate 
obtained in Sec. 4.1. We determine the magnetic mo- 
ment of each magnet by using the susceptometer to 
measure the force between all possible pairs of magnets, 
aligned coaxially and placed a known distance apart. 

We first set the span of the susceptometer so that Zo 
is approximately 15 mm. A precise value of Zq is deter- 
mined from a measurement of F^ using the Alacrite 
standard whose calibration was described in Sec. 4.2.1. 
We then add gauge blocks with a height Zb = 70 mm so 
that Zo is nominally 85 mm. 

At this setting, we measure the force between mag- 
nets A and B. Magnet B is placed on the span so that the 
resulting interaction force is attractive. Its orientation is 
then reversed so that the reaction force is repulsive. The 
same sequence of measurements is repeated for C. En- 
suring that the two magnets are sufficiently coaxial is 
relatively straightforward because the coaxial condition 
gives a maximum change in the balance reading. 
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The average magnitudes of the attractive and repul- 
sive forces give us /% between magnets A, B and A, C. 
We then place magnet B on the pedestal and repeat the 
measurements by placing A and C on the span in each 
orientation to obtain /% between B, C and B, A. The last 
measurement serves as a check on the reproducibility of 
the results. As may be inferred from Table 2, we found 
the balance readings were nominally 8 300 |jug. The 
mean standard deviation of four repeated balance read- 
ings was 15 |jug. 

The force measured between magnets A and B is 



(Fb)AB = Fa 



4lT 



morris 
{Zo + LI2f 



where LI2 is 2.5 mm, half the height of the magnet, and 
we have added subscripts to distinguish the moments of 
magnets A and B. The equation is the well known rela- 
tion for the interaction force between two coaxial 
dipoles and is a special case of Eqs. (5a) and (6b) when 
Zo is large compared with the sample dimensions and x 
is negligible. Thus 



2 f AB ■ f AC ATTiZo + LITf 

m-A = — 7^ 



6^io 



(7) 



The indices A, B, and C in Eq. (7) can be permuted 
in obvious fashion to obtain the moments of magnets B 
and C. 

The results are given in Table 2. We note that, on 
average, the three magnetic moments agree well with 
the estimate of 0.0898 A-m^ obtained in Sec. 4.1. 



Table 2. Calculation of magnetic moments m of magnets A, B and 
C based on the force Fi, measured between pairs of magnets aligned 
at a spacing of 84.90 mm 



Magnets 



fb/(Ji-N° (measured) 



m/(A-m ) (calculated) 



A, B 


81.80 


A, C 


83.35 


B, C 


83.50 


A, B 


81.86 


A 




B 




C 





0.08910 
0.08937 
0.09089 



' +1 [lN corresponds to a balance reading of -101.9 [Lg. 



With the value of m in hand, we can now compute x's 
for the Alacrite standard. The result is 



X\ = 0.001 348 

with MC = 0.000 012 (Type B standard uncertainties 
neglected). 

4.3 Calibration Checks 

Once we have a standard with a known susceptibility, 
any height Zo is determined by a measurement of f a(s), 
where (s) denotes use of the standard. The susceptibility 
X' of an object of interest is then determined through a 
measurement of the force F^ due to the unkown object 
placed at the same Zq. In essence, the following relation 
is used: 






fa /a(s) 

Fa(s) /a 



(8) 



where /a(s) is appropriate to the standard and /a to the 
unknown [Eqs. (5b) and (6a)]. In Sec. 3.2 we referred to 
the special case where the standard and unknown have 
the same dimensions so that {/a(s)//a} = 1 in general. 
This relation also holds for "large" samples as defined in 
Sec. 3.1. The ratio of susceptibilities is manifestly insen- 
sitive to details of the theory for such cases. It is only 
when the term { /a(s)//a} depends sensitively on Zo that 
we must rely on the quantitative validity of the theory 
developed above. 

Two measurements described in detail in Ref. [1] 
provide checks on the method of calibration used in the 
previous section and help in the assessment of overall 
type B standard uncertainties. First, it was shown that 
the susceptibility measured for a number of well-char- 
acterized materials is generally in good agreement with 
handbook values. By itself, this is not a sufficient check 
of the apparatus because the susceptibilities of the stan- 
dard materials used are relatively small in magnitude 
compared with the susceptibilities of stainless steels and 
other alloys of interest. Thus the standard materials can 
only be used to verify operation at relatively high fields, 
corresponding to values of Zo in the range of 1 mm to 
15 mm. In this range, the approximation of the magnet 
as a dipole is expected to bias results for samples that are 
not "large." We expect good results for large samples, 
however, and this serves to check our inferred value of 

X\- 

The second test is that the ratio of the measured 
susceptibilities of two different LIH alloys having sig- 
nificantly different dimensions and susceptibility is 
found to be independent of Zo through a range of settings 
where {/a(s)//a} changes by 50 %. These tests suggest 
that, for the examples given below, the type B relative 
standard uncertainty is no greater than 3 % in the mea- 
sured susceptibility for LIH samples. 
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5. Typical Use of Susceptometer 

We find it convenient to use the automatic-zero fea- 
ture of the balance. The reading displayed is conse- 
quently insensitive to small, slowly changing forces. We 
therefore take the following precautions during the mea- 
surements. We do not collect data at settings where the 
introduction of the sample changes the balance reading 
by less than 10 counts in the last displayed digit. To 
change the position of the sample, we remove it entirely, 
let the balance return to zero and then place the sample 
in its new position. We repeat a measurement if the 
balance zero changes significantly during a sequence of 
readings. 

For routine measurements, we begin with Zo at no less 
than 25 mm and only decrease Zo if the sample suscep- 
tibility is too small to produce a reasonable signal. Set- 
ting Zo to 25 mm corresponds to a maximum field 
strength of H^^ ~ 900 A-m"' {fjbvnH^^^ ~ 1 mT). 

Appendices B and C summarize susceptibility mea- 
surements for typical 1 kg mass standards of Class E2. 
The first (Appendix B) has the shape of a right-circular 
cylinder with rounded edges and the second (Appendix 
C) has the external dimensions of a Class M standard 
[6]. 

5.1 Cylindrical Mass Standard 

In the example shown in Appendix B, the sample was 
placed coaxial with the magnet at a nominal spacing of 
30 mm. Measurements were made using both vertical 
orientations of the magnet and both vertical orientations 
of the sample. Data obtained using the Alacrite standard 
were used to determine the average Zq, 30.55 mm. Equa- 
tions (5a) and (6b) predict that the Alacrite standard 
should produce a value of F\,l g = +5 [xg. We observe 
+6 [xg, which is a useful check that the device is operat- 
ing properly. 

As shown in Appendix B, the balance readings ob- 
tained when the bottom of the sample was nearer the 
magnet differ from those when the usual orientation of 
the sample is reversed. This could be a real difference 
due to sample inhomogeneity. However, the fact that Fa 
is essentially identical for each orientation of the sample 
suggests that permanent magnetization is the more prob- 
able explanation. 

We emphasize that our calculation of M^ assumes a 
model of uniform magnetization. Evidence that this 
model is unrealistic is that the fitted value of M^ for the 
sample in its normal orientation is 0.2 A-m"' but 0.07 
A-m"' when the sample is reversed. 

Based on the limit derived in Sec. 2, we would not 
expect the observed magnetization to be of concern. To 
see if the sample magnetization has any practical 



consequences, we measured the mass of the secondary 
standard using a magnetically servocontrolled balance 
(Mettler-Toledo HK 1000 MC). The results obtained 
were the same when the standard was upside down on 
the balance. This means that, although the sample has a 
detectable permanent magnetization, routine weighings 
in our laboratory should be unaffected. 

5.2 OIML-Shaped Mass Standard 

Mass standards are not generally cylinders, although 
they do possess cylindrical symmetry. A typical shape 
is the so-called OIML design [6], shown in Appendix C. 
It differs from a simple cylinder primarily in that it has 
a lifting knob on top and a recessed base. 

Appendix C presents calculations of ;^'' for a 1 kg 
standard in its usual orientation carried out at a setting 
of Zq typical for routine measurements. Both orienta- 
tions of the magnet are used in order to obtain the value 
of Fa. The same data are used to infer an apparent value 
of Mj from the calculation of F^, see Eq. (5). 

To demonstrate the generality of the theory, we also 
made measurements with a similar mass standard 
placed on its side so that its body was centered directly 
over the magnet. Equation (5b) was solved numerically, 
with limits of integration corresponding to the cylindri- 
cal body of the standard. As shown in Table 6 (Appendix 
C), the presence of the knob has a relatively small effect 
on the calculation. Experimental results were un- 
changed when the sample was rotated about its axis of 
symmetry. 



6. Behavior of Nonlinear and/or Inhomo- 
geneous Samples 

The previous section described typical measurements 
for alloys used in the manufacture of good secondary 
mass standards. We now give examples of results ob- 
tained from test objects that are nonlinear and/or inho- 
mogeneous. 

6.1 Nonlinear Alloy, Initially Unmagnetized 

Pure brass is diamagnetic. We find that much of the 
industrial brass supplied to our workshops is, on the 
contrary, contaminated with magnetic impurities to the 
extent that the measured susceptibility is frequently 
greater than 0.01 in field strengths ranging up to about 
1 kA-m"'. Typically, the susceptibility depnds on maxi- 
mum field strength. 

If one attempts to make the measurement at a maxi- 
mum field strength of 10 kA-m"' or more, the balance 
reading does not stabilize in the usual settling time but 
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continues to increase in magnitude for many minutes. 
This is an indication that a portion of the sample is 
becoming permanently magnetized, as can be verified 
by carrying out further measurements with a maximum 
field an order of magnitude less strong. Exposure of 
poor-quality samples to high fields must, therefore, be 
avoided. This is the reason we begin all routine testing 
at low fields. For the same reason, we take great care 
that test objects do not accidentally come in contact with 
the permanent magnet. 

6.2 Nonlinear Alloy, Initially Magnetized 

It is not unusual to encounter samples of industrial- 
grade copper alloys or of type 304 stainless steel that 
have accidentally been magnetized at some point in their 
history. Once magnetized, stainless steels are difficult to 
demagnetize, as is well known [5, 10]. 

We found one disc of copper-beryllium alloy (mass 
fraction of beryllium nominally 2 %) that was magne- 
tized to such an extent that the term in M^ in Eq. (5a) 
dominated the measurements. Results for this sample 
are given in Table 3. It may be seen that the model of 
uniform magnetization fits reasonably well in this case, 
with M^ about 25 A-m"'. Commercial alloys of 
Cu-2 %Be contain small but significant amounts of 
ferromagnetic elements. It is, therefore, the heat treat- 
ment given to the alloy that renders it more or less 
magnetic [11]. 

We then passed the sample through a demagnetizing 
coil energized at the main frequency (50 Hz). Results 
obtained after this operation are also given in Table 5. 
We see that, although still present, the permanent 
magnetization is weaker and evidently not as uniform. 
Data obtained at the two nearest settings of Zo now 



reveal a measurable F^, indicating a susceptibility of 
order +0.000 35. This value is consistent with the range 
of susceptibilities common to commercial alloys of 
Cu-2 %Be [11]. 

6.3 Good-Quality Alloy, Surface Inhomogeneity 

One of the well-characterized materials that we tested 
was a disc of oxygen-free copper with a mole fraction 
for iron impurity of only 2X10"^. This alloy has a hand- 
book susceptibility of -9.6X10'^ [11]. The disc-shaped 
sample supplied by the manufacturer had rough-sawn 
faces that we machined flat before the measurements. 

The first susceptometer reading appeared to show 
that the sample was paramagnetic. After it was given a 
light etch in a copper-cleaning solution, further mea- 
surements showed the sample to have the expected sus- 
ceptibility. Thus the initial problem was a surface effect, 
presumably due to iron impurities transferred from the 
cutting tool. In order to have sufficient signal, these 
measurements were made at Zq = 10 mm. 

Hard-working of stainless steel surfaces, as may 
happen during the fabrication of mass standards, can 
cause the susceptibility to increase [4]. Should this 
occur, the fabricated standard becomes inhomogeneous 
and use of Eq. (5b) will make it appear that the suscep- 
tibility of the sample is decreasing with increasing Zg. 
This may be the explanation of results reported in Ref . 
[1] for a cylinder of the stainless steel Immaculate V 
where the susceptibility inferred from Eq. (5) decreased 
by about 7 % as Zq was increased from 9.3 mm to 24 
mm. The susceptibility did not decrease significantly 
when Zo was further increased to 29 mm. The sample 
had no detectable permanent magnetization. 



Table 3. Effect of "degaussing" on a magnetized disc of commercial-grade Cu— 2 %Be. The influence of the 
Earth's field on the calculation of Mj may be neglected 





Before degaussing 




After degaussing 






Nominal 


Measured 




Measured 




Measured 






Zo/mm 


Fb/hLN" 


MJ{A-m-') 


FJ[i.N 


A-ZIO^ 


FJ[x.N 


M, 


/(A-m"') 


15 






0.96 


3.52 


7.51 




13.2 


20 


12.68 


25.5 


0.27 


3.40 


5.80 




12.1 


25 


10.06 


25.2 












30 


7.72 


24.9 






3.14 




10.7 


40 


4.395 


24.4 












50 


2.531 


24.0 






0.97 




9.7 


70 


0.963 


24.0 






0.37 




9.7 


90 


0.404 


22.8 













"+1 [lN corresponds to a balance reading of —101.9 \xg. 
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We cannot take quantitative account of an enhanced 
surface susceptibility without prior knowledge of either 
the thickness of the enhanced layer or its susceptibility. 

7. Magnetic Fields in a Balance 

It is possible to measure both the vertical magnetic 
field intensity and its gradient in the weighing chamber 
of the MT5 balance used for these measurements. Using 
a conventional magnetometer, we find the vertical field 
within the weighing chamber (with the cylindrical mag- 
net removed) is about 10 % greater than that of the Earth 
alone. 

The gradient of this field, measured by raising the 
probe a few centimeters, was inferred to be of order 
40 A-m"^. An independent measure of the gradient is 
obtained by observing the change in balance reading 
when the magnet is reversed with no sample present. 
The measured difference is 1.5 mg implying a gradient 
of 65 A-m"^ (note that the difference in balance reading 
upon reversal of the magnet reported in Ref. [1] was 
0.9 mg; this was obtained from a different unit of the 
same model balance). 

Thus, based on Eqs. (2) and (3), if the balance is used 
to weigh 5 g of stainless steel, the susceptibility of the 
material should be less than 5 and any uniform magne- 
tization in the vertical direction should be less than 
200 A-m"' in order to keep magnetic effects below the 
balance resolution of 1 |jug. If the same average field 
strength and gradient were present in a 1 kg mass com- 
parator, stainless steel standards would require 
X < 0.025 and IM^I < 1 A-m"' in order to keep magnetic 
effects below 1 [xg. 

8. Conclusion 



samples. An acceptable upper limit for uniform magne- 
tization has been derived in the spirit of a previous 
discussion of an acceptable upper limit for the suscepti- 
bility 



Appendix A. 
paratus 



Construction of the Ap- 



The apparatus is shown schematically in Fig. 1 and 
the most important construction details are given below. 

9.1 Balance 

The balance used (Mettler-Toledo AG, model MT5) 
has a 5 g capacity and a resolution of 1 |jug. The weigh- 
ing chamber is a glass tube with a removable glass top. 
We replace the original top with one made of an alu- 
minum alloy, as shown in Fig. 6. The alloy we use, 
designated AU4G in France, corresponds closely to 
ASTM 2017 and is the aluminum alloy commonly 
stocked by our workshop. 




Top 



We have discussed how magnetic forces may lead to 
errors in mass measurements. The force equations that 
describe the unwanted effects also describe the opera- 
tion of the BIPM susceptometer. We have presented 
design details of this device, which is suitable for check- 
ing the susceptibility of 1 kg mass standards of stainless 
steel and similar nonmagnetic or weakly magnetic al- 
loys. 

In addition, we have demonstrated that the device may 
be calibrated by using a set of gauge blocks, two addi- 
tional magnets, and a LIH sample of unknown suscepti- 
bility. The calibration was verified by measurements, in 
relatively high field strengths, of materials having well- 
characterized susceptibilities. 

The device can also detect a permanent magnetiza- 
tion. An apparent magnetization can be calculated and 
has a simple interpretation for uniformly magnetized 



80 



74 



Front 

Fig. 6. Replacement cover for the balance, made by gluing an annu- 
lar piece of aluminum alloy (AU4G) to the bottom of a thin disc of the 
same material. Dimensions are in millimeters. 



9.2 Magnet 

The cylindrical magnets used are made of 
neodymium-iron-boron (Vacuumschmelz GmbH, type 
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370 HR). The dimensions of each magnet are height = 
2.5 mm; diameter = 5 mm. The axis of magnetization 
coincides with the geometric axis of the cylinder. The 
density of the NdFeB alloy is 7400 kg-m"^. Two such 
magnets are combined to produce a cylinder of height 
equal to diameter. In all equations given in the text, the 
origin of the coordinate system is the geometric center 
of the magnet. 

When in doubt, the north pole of the magnet can 
easily be located by placing the magnet horizontally on 
a smooth, nonmagnetic surface. It then acts like a com- 
pass, spontaneously aligning itself along the Earth's 
horizontal field. 

9.3 Pedestal 

The magnet sits on a tubular pedestal of AU4G that 
is itself centered on the balance pan. Pedestal dimen- 
sions are shown in Fig. 7. Holes (not shown) drilled 
through the wall of the tube ensure that the total mass of 
the magnet and pedestal is well within the capacity of 
the balance. 

An aluminum shim can be used to raise the magnet 
closer to the balance cover so that Zo is reduced to a 
minimum, thereby achieving maximum sensitivity. The 
shim is not used in routine work. The total mass of the 
magnets, pedestal and shim is about 4.2 g. 

9.4 Bridge 




Front 



Samples are centered on a bridge, also made of 
AU4G, that straddles the weighing chamber of the bal- 
ance. Details are shown in Fig. 8. The important fea- 
tures of the design are that the span can be made level 
with respect to the balance pan and that the thickness of 
the span is the minimum consistent with adequate me- 
chanical rigidity. If the span sags in the middle, samples 
of large diameter will be farther from the magnet than 
samples of small diameter. For critical measurements, a 
diamagnetic shim of small diameter can be placed at the 
center of the span thereby overcoming the problem of 
sag. 



9.5 Gauge Blocks 

We have found it convenient to place the bridge on 
nonmagnetic gauge blocks so that it may be raised or 
lowered in precise increments of 5 mm. We manufacture 
our blocks either of AU4G or of selected brass. The 
selection procedure consists of placing the brass stock in 
the vicinity where the eventual gauge block will be 
placed and verifying that there is no change in the bal- 
ance reading with the magnet installed. 



Fig. 7. Pedestal used to support the magnet. Not shown are 48 holes 
(eight columns of six), 5 mm in diameter, that are drilled into the tube 
in order to reduce its mass. The material used for the pedestal is 
AU4G; dimensions are in millimeters. 



9.6 UMTS Balance 

A limited set of measurements was made using a 
model UMT5 balance (Mettler-Toledo AG), which has 
the same capacity and overall dimensions as the MT5 
but has a resolution of 0. 1 (xg. The pan supplied with the 
UMT5 has a diameter of only 16 mm and, in order to 
use the pedestal shown in Fig. 7, the original pan was 
replaced by the pan from an MT5. Results obtained 
were satisfactory and indicate that the increased resolu- 
tion allows the determination of the susceptibilities of 
samples of correspondingly smaller volumes. 
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Fig. 8. Bridge used to support the test samples. Not shown are 
screws that are used to fasten the horizontal span to the two vertical 
supports. Also not shown are fiducial marks lightly inscribed on the 
top surface of the span and used as an aid in centering samples. The 
material used is AU4G except for the screws, which are brass. Dimen- 
sions are in millimeters. 



10. Appendix B. Typical Calculations for 
a Cylindrical Mass Standard Coaxial 
with the Magnet 

The following describes how we routinely determine 
the susceptibility and effective permanent magnetiza- 
tion of a cylinder oriented coaxially above the magnet 
(Fig. 1). We find it convenient to work with spreadsheet 
software that calculates: Fa, from Eqs. (4), (5) and (6) as 
well as Fb, from Eqs. (5) and (6). The value of F^ is 
displayed as two separate components: from the Earth's 
vertical magnetic field Hez and a possible permanent 
magnetization M^ coaxial with the sample. As we 
discuss below and in Sec. 1 1 , it is also useful to display 
/a calculated from Eq. (5a). This is a number between 1 
and that takes account of the finite dimensions of the 
sample relative to a semi-infinite slab. Spreadsheet 



variables are the radius c and length t of the cylindrical 
sample, Zo, x^ ^nd M^ . 

To find Zo, we start with the measurement of Fa made 
with the Alacrite standard. With c, t and x\ appropriate 
to this sample, we call upon a spreadsheet utility to 
compute the value of Zq, through iterations of Eqs. (5a) 
and (6a), forcing the calculated value of Fa to agree with 
that found by experiment. We check that Fb with M^ = 
is consistent with the computed result. 

We then enter c and t appropriate to the sample under 
test and, calling once more on the spreadsheet utility, 
adjust X until the measured value of Fa is obtained. If 
necessary, we also adjust the value of M^ until the calcu- 
lated force Fb agrees with experiment. A separate calcu- 
lation for X and M^ is carried out for the sample in its 
normal orientation and turned upside down. 

Typically, the sample under test has edges that are 
rounded with a radius of about 2 mm. To check whether 
the actual shape produces the same signal as a cylinder 
with sharp edges, we make supplementary calculations 
shown in Table 4. 

The simplest calculation, which is usually sufficient, 
is to set limits on the susceptibility. That is, the suscep- 
tibility must be less than that computed assuming a 
perfect cylinder of radius c and length t . The suscepti- 
bility must be greater than that computed for a perfect 
cylinder of radius (c-2 mm) and length t. The results of 
these calculations are shown in Table 4. It is sufficient to 
calculate /a since x is then found directly from Eq. (8). 

We can also calculate the susceptibility more pre- 
cisely. By superposition, we first determine I^ due to the 
smallest annular volume that contains the lower rounded 
edge. The contribution is, in fact, small compared with 
the total /a. We may be satisfied with this approximation; 
or we may note that all parts of the annulus are approx- 
imately equidistant from the magnet so that the signal 
from any arbitrarily chosen region of the annulus is 
roughly proportional to the volume of the region, as may 
be inferred from Eq. (5b). In particular, the material cut 
away from the annulus to make the rounded edge has a 
volume that is (l-7r/4) times the total annular volume. 
Thus we have underestimated x by about 0. 1 %, an error 
that is well below the uncertainty of the final result. The 
correction for the upper rounded edge is even smaller, so 
we neglect it. 

11. Appendix C. Typical Calculations 
for an OIML-Shaped Mass Standard 

The mass standard is shown in Fig. 9. The lifting knob 
and the recessed base are complications that are dealt 
with in detail in this section. Two orientations of the 
sample are now considered. 
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Table 4, Computation of x and Mj for a 1 kg stainless steel standard of cylindrical shape. The 
height and diameter of the cylinder are both 54 mm and the edges are chamfered with a radius 
of 2 mm. Measurements were made with the cylinder coaxial to the magnet at a spacing of 
Zo = 30.57 mm 



Experimental data 


Sample orientation 


F./H.N° 


Fb/lii-N 


Normal 
Reversed 


0.902 
0.904 


0.489 
0.199 


Calculation of x 


Sample dimensions 


/. 


X 



Height = diameter = 54 mm 


0.773 




0.006 85 


Height = 54 mm; diameter = 50 mm 


0.733 




0.007 21 


Height = diameter = 54 mm 


0.773 






Height = 2 mm; diameter = 54 mm 


-0.199 






Height = 2 mm; diameter = 50 mm 


+0.192 






Correction for lower edge 


-0.007 (1-17/4) = 


-0.002 




Total 


0.771 




0.006 86 



Calculation of Mj 


Sample orientation 


Measured 
Fb/fJi-N 


Theoretical 
effect of H,^/jjlN 


Adjusted 
Mj/(A ■ m-') 


Normal 
Reversed 


0.489 
0.199 


0.275 
0.275 


0.20 
0.07 



"+1 (jlN corresponds to a balance reading of —101.9 (JLg. 



11.1 Sample Coaxial With the Magnet 

The usual experimental orientation we use when de- 
termining magnetic susceptibility is coaxial with the 
magnet. We first measure Zo as described in Appendix 
B. It is then possible to put simple limits on the suscep- 
tibility: it must be less than that based on a sum and 
difference of "outer" cylinders that contain the entire 
volume of the standard; it must be greater than that 
based on a sum and difference of "inner" cylinders that 
are contained within the mass standard (see Fig. 9). The 
results are given in Table 5. Note that the net contribu- 
tion of the knob is virtually negligible. 

The same value of Zo was used in calculating each 
component of l^. In general, any component l\ = h{Z'o) 
calculated using a different distance Z o must be normal- 
ized by a factor (Zq/Z o)** before combining it with other 
components calculated using Zq. The reason for this is 
easy to see: Consider an example where there are only 
two components in the summation, I^ and /'a , computed 
at distances Zo and Z'o respectively. From Eq. (5a) and 
the principle of superposition, the total measured force 
must be described by the relation 



"T^ 




rr 








-T-r 


) ( 
















1 




2 




3 




4 



"OUTER" CYLINDERS 




^ 



( 




) 















rr 



2 3 

INNER " CYLINDERS 



VX' ■Fmax(Zo)] • h+[X' ■Fmax(Z'o)] ' /' 



(9) 



Fig. 9. Outer and inner cylinders used to approximate an OIML- 
shaped 1 kg mass standard by the method of superposition. Dimen- 
sions of the cylinders are given in Table 5. 
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Table 5. Computation of x and M, for an OIML-shaped 1 kg stainless steel standard. Details of the shape are 
given in Fig. 9. Measurements were made with the cylinder coaxial to the magnet at a spacing of Zo = 24.70 mm 









Measured values Fa =1-034 
Fb = 0.169 


p,N 
H,N 








Calculation of x 








c/mm 


f/mm 






/„ 


X 


Outer cylinders 
1 

2 
3 
4 






24.0 
16.5 
21.5 
21.5 


58.5 

1.0 

80.5 

58.5 

Total 






0.8196 
0.1116 
0.7690 
-0.7681 
0.7089 


0.00364 


Inner cylinders 
1 

2 
3 
4 






24.0 
17.5 
13.5 
13.5 


58.5 

1.0 

80.5 

58.5 

Total 






0.8196 
-0.1164 

0.4979 
-0.4975 

0.7036 


0.00367 


Calculation of Mj 


Sample dimensions 








Measured 
Fb/(xN 


Theoretical 
effect of //e,/p.N 


Adjusted 
Mj/(A-m-') 


OIML 

Cylinder, diam. = 48 mm 
correction for cylindrical 
Correction for bevel 


ht. = 59 mm 
recess 

Total 


0.169 




0.199 
-0.012 
-0.001 

0.186 




S0.07 



' +1 jjlN corresponds to a balance reading of -101.9 \xg. 



By simple algebra, 



/a + 



Zo 



/'a 



(10) 



Equation (10) can easily be generalized to an arbitrary 
number of components. 

Once the total /a has been found, we can determine x 
from Eqs. (4) and (5a). Since F^, is a relatively small 
number, we may approximate the sample as a cylinder 
resting on an annulus for this calculation. A more de- 
tailed calculation, shown in Table 5, verifies this asser- 
tion. 

11.2 Axis of Sample Perpendicular to Axis of 
Magnet 



coordinates. The computation can be done using pack- 
aged software such as Mathcad: 

Zon = 1 



e = 10 ' 

/a = - 



3tt 



d 
.3Zo„ 



/-Zo„+2c„ ^ V c2_(Zo„ + c„-z)2 + e 
•'Zn„ •'o 






2 , ^,2 , „2\A dxdydz 



It is possible to measure the susceptibility of a sample 
laid on its side, its body centered about the magnet axis. 
Our approach is, again, to calculate I^ and then solve for 
X using Eqs. (4) and (5a). In this case, however, /a is 
determined from Eq. (5b) by numerical differentiation 
of the volume integral, after first converting to Cartesian 



All integration limits are dimensionless, having first 
been normalized to Zo = 25.39 mm. The small quantity 
e is added to ensure that the upper limit of the second 
integral is real. Integration is over the dummy variables 
X , y and z representing the two horizontal axes and the 
vertical axis of a cartesian coordinate system where the 
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Table 6. 


Calculations for horizontal 


orientations 


of the mass 


standard 








c/mm 


f /mm 


h 


X 


24 


59 


0.661 


0.003 68 


24 


21 + 59 + 21 


0.684 


0.003 55 


24 


21+59 


0.672 


0.003 62 



X-axis is parallel to the sample axis. Computing time can 
be shortened considerably by analytical integration over 
the X -variable, a result that is found in standard refer- 
ences on integral calculus. The results of the calculation 
are shown in Table 6. The measured value of the attrac- 
tive force was F^ = 0.871 (xN. The first calculation ex- 
cludes the knob in the calculation of I^ and therefore 
overestimates x- The second assumes there is a knob on 
each end of the body of the standard and that each knob 
is a simple extension of the body by 20 mm; this under- 
estimates X- The final calculation (the average of the 
previous two) assumes a single knob but models it again 
as a simple extension of the body. This again overesti- 
mates X t)ut gives a tighter bound. Therefore 
0.003 68 >;f> 0.003 62. 

If the sample were permanently magnetized perpen- 
dicular to its axis of symmetry, the measurements of F^ 
would depend on the azimuth of the sample and this can 
easily be checked. It is, of course, possible to estimate 
by numerical methods the expected value of Fy, due 
simply to the Earth's field. 

The measurements summarized in Tables 5 and 6 
were made on different 1 kg standards fabricated from 
the same alloy. 
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[J. Res. Natl. Inst. Stand. Technol. Volume 100, Number 3, May- June 1995, p. 209] 



Page 211, above Eq. (3). 

Replace "the magnetic field within the balance 
chamber is vertical with the form..." by "the vertical 
magnetic field strength within the balance chamber has 
the form...". 

Page217, Eq. (8). 
Replace "/a" by "/,". 

Page 218, two paragraphs above Sec. 5.1. 
Replace "//m//^^;' by "//(//„^;'. 

Page 223, first paragraph of Sec. 11.1. 

Replace ": it must be less than that based on a sum...; 

it must be greater than that based on a sum..." by 

": it must be greater than that based on a sum...; it must 

be less than that based on a sum...". 



Page 224, Table 5, column 4, row 2. 
Replace "0.1116" by "-0.1 11 6". 



Page 224, Eq. (10) should be: 



(z ^ 
z' 

V V 



r 



Page 225, two paragraphs above the Acknowledgment. 
Replace "This again overestimates % " by "This again 
underestimates x'\ 
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